Theory of Tunneling Spectroscopy in a Mni2 Single-Electron Transistor by 

Density-Functional Theory Methods 
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We consider tunneling transport through a Mni2 molecular magnet using spin density functional 
theory. A tractable methodology for constructing many-body wavefunctions from Kohn-Sham or- 
bitals allows for the determination of spin-dependent matrix elements for use in transport calcula- 
tions. The tunneling conductance at finite bias is characterized by peaks representing transitions 
between spin multiplets, separated by an energy on the order of the magnetic anisotropy. The energy 
splitting of the spin multiplets and the spatial part of their many-body wave functions, describing 
the orbital degrees of freedom of the excess charge, strongly affect the electronic transport, and can 
lead to negative differential conductance. 
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There is a growing interest in exploring the rich physics 
and spintronics functionality of molecular single-electron 
transistors (SETs) consisting of a few magnetic molecules 
weakly coupled to nano-gapped electrodes Recently 
two groups 0, Q have carried out single-electron tunnel- 
ing experiments on individual magnetic molecules based 
on Mni20i2 (henceforth Mni2) with organic ligands. 
Mni2 is the most studied and perhaps the most remark- 
able molecular magnet In its crystal phase, Mni2 
is characterized by a long spin relaxation time due to 
its large uniaxial magnetic-anisotropy energy. Further- 
more, at low temperatures, quantum effects in the re- 
laxation properties are clearly discernible and have 
been attributed to quantum tunneling of the molecule 
collective magnetization How these properties are 
revealed in electronic quantum transport is a question 
of CTeat significance for the field of molecular spintron- 
ics d^. Indeed, the SET experiments 0, Q show signa- 
tures of the molecule magnetic state and its low-energy 
collective spin excitations. The theoretical models pro- 
posed so far , are typically based on effective 
giant-spin Hamiltonians with large uniaxial anisotropy 
barriers. This approach has two drawbacks [16]. First, 
the effective spin Hamiltonian for the charged states 
(anion and cation) of Mni2 needed to describe sequen- 
tial tunneling transport, is not known. Scaling of the 
global anisotropy parameter to account empirically for 
changes in the number of electrons forming a macro- 
spin is fraught with uncertainty [l7| . Second, the or- 
bital degrees of freedom are not included in the giant- 
spin Hamiltonian formalism. The orbital effects due to 
changes in electron population on the Mni2 molecule 
modify the symmetry and magnitude of the magnetic- 
anisotropy Hamiltonian and can even change the spin 
ordering [l8|. 

In this Letter we provide a microscopic many-body de- 
scription of the ground state (GS) and low-lying spin 
excitations of both neutral and charged states of a Mni2 
molecular magnet. Our approach is based on spin density 



functional theory (SDFT), which has been very successful 
in describing the spin-orbit-induced magnetic anisotropy 
barrier in Mni2 and other molecular mag nets [ii-ilj. We 
find that when a delocalized electron is added to (sub- 
tracted from) the molecule, the GS spin of the molecule 
increases (decreases) by 1/2. For both charged states, 
the GS magnetic-anisotropy energy is larger than for 
the neutral Mni2. We then incorporate this information 
into a quantum master equation for electronic transport 
in the sequential tunneling approximation, which is ap- 
propriate for the experimental Coulomb blockade (CB) 
regime. The approximate many-particle eigenstates lead 
to a tunneling conductance that exhibits fine structure 
on the order of the anisotropy energy and, under cer- 
tain circumstances, to strong negative differential con- 
ductance (NDC) . Comparison with the giant-spin model 
shows that spatial selection rules play a crucial role in 
determining which spin excitations contribute the most 
to the tunneling conductance. 

We need to know the many-electron wave functions, 
representing low-energy spin excitations, as a function of 
the excess charge (Q), spin ordering (M = {iJiv}u=i,i2)i 
applied electric E and magnetic B fields, and the pa- 
rameters 9, (j) describing the quantization axes. We refer 
to the collection of all possible variables as the "order- 
parameter vector" (OPV), p = (Q, M, E, B, 6*, 0) to la- 
bel the states. Given a specific OPV, we first construct 
a set of Kohn-Sham (KS) single-particle states $fe(p), 
by diagonalizing a KS single-particle Hamiltonian 77 (p) 
that depends upon this OPV. Since some of these effects 
(Q, M) are clearly large and some (6*, (j)) are generally 
small, there is flexibility as to which of these terms must 
be accounted for self-consistently. Specifically, 



i/(p) = i/(Q,M,E,B,0,, 
= i/o(DFT,Q,M)-l 
+B- (L-K2S), 



E 



(1) 



contains a spin-polarized term i/o(DFT, Q, M), which is 
treated self-consistently for the cation, neutral and an- 
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ionic states {Q = +1,0 and —1); Vl s(^, 0) represents 
the spin-orbit interaction. We neglect the last two terms 
representing the coupling to external fields. 

The spin-ordering M corresponds to that obtained 
from the local moments of the 12 Mn atoms (fj,^) in 
the classical ferrimagnetic state of the neutral molecule 
22 1 . The spin-orbit operator is treated exactly [31 ~ 
albeit non-self-consistently - in the basis of the eigen- 
states of i7o(DFT, Q, M). Diagonalizing the above 
Hamiltonian with the constraint that the expectation 
value of the total spin ((S)) is quantized along the 
axis determined by 6,(j) results in a set of single- 
particle, noncollinear spin orbitals, 0fc(p), expressed 
as (j)k = (l)l{r)x+{0,(t>) + <t>l{r)x-{d,4)). Here, the 
angle-dependent spinors x±(^i0) spin-1/2 coherent 
states specified by the quantization axis, x±{^i4>) = 
cos(6'/2)|±) ± sin(6'/2)e±''^|T)- 

We now construct approximate many-body functions 
for the ground and excited electronic states as single 
Slater determinants (SDs) of the spin orbitals 4>k{p)' 



\p;ki,k2,. . .kNg) = |(?!)fci(p)(/)fc2(p)...(/)fc„^ (p)). 



(2) 



The above states are generally not eigenstates of either 
or 5*2. However, a state with 1(8^)1 = So{So + 1), 
especially when constructed from a closed shell of spa- 
tial states, is expected to be the primary contributor to 
an 5 = 5*0 eigenstate. While the variables {6, (/>) gener- 
ate a continuous overdetermined set of SDs, a judicious 
choice of 250 + 1 values of 9 and (f> can lead to a nearly 
orthogonal set of normalized linearly independent many- 
particle SDs, with (Sz) taking on integer or half integer 
values akin to the standard M = — 5*0, — -|- 1, . . . , 
quantum numbers. Choosing 25-1-1 values of given by 
Scos{9]\f) = M leads to integer or half integer {Sz) re- 
gardless of the choice of 0; however, choosing 4>m = M(j)o 
leads to destructive interference in the off-diagonal el- 
ements of these states and aids in producing approxi- 
mate Sz eigenstates. For the case of S* = 10, we find 
that choosing 0o — 34° leads to the smallest off-diagonal 
overlaps between approximate eigenstates. We call the 
2^0 -t- 1 many-electron states constructed with this pro- 
cedure a spin multiplet. Besides the GS spin multiplet, 
the anion and cation have a few low-lying excited spin 
multiplets. These come about because the HOMO level 
of the charged molecule is quasidegenerate with a many- 
fold of LUMO levels [18[ . Using Eq. (0) , we can construct 
several SD excited-states close in energy to the GS, all 
having the same spin, Sq- The relevant spin- multiplets 
for the molecule are shown in Fig. [TJa). Typically, the 
level spacing within each spin multiplet is of the order 
of 0.1-1 meV, while different multiplets are separated by 
energies of the order of 10 meV. Note that the energies 
of a given spin multiplet are not exactly invariant un- 
der M —M, since the choice of the angle (j)M is in- 
commensurate with the nonperfectly uniaxial symmetry 
of Mni2. The breaking of the level degeneracy for the 



(Af , —M) pair is of the order of the transverse anisotropy 
terms coming from 4th order spin-orbit contributions, 
and therefore very small (2^. For a later comparison 
with the giant-spin approach, we disregard the small de- 
viations from uniaxial symmetry and consider 



TABLE I. The GS properties from DFT: spin, energy and 
magnetic-anisotropy energy as a function of charge. 



State 


Q 


Spin 


Energy (eV) 


MAE (K) 


MAE (meV) 


Anion 


-1 


21/2 


-3.08 


137 


11.8 


Neutral 





20/2 


0.00 


55 


4.7 


Cation 


1 


19/2 


6.16 


69 


5.9 



the spin Hamiltonian Hsq = J2i,n ^SQ,i,n[SQi]"' , given 
in terms of spin variables Sqa for each spin- multiplet 
i of charge state Q; DsQ,i^n are anisotropy constants 
which we extract by fitting the corresponding SDFT 
energy spectra. Table |T] shows the total spin, GS en- 
ergy and magnetic-anisotropy energy of the lowest spin 
multiplet calculated for the different charge states of 
Mni2. Adding an electron to the molecule increases the 
molecule spin and decreases the energy. Furthermore, 
the anisotropy increases significantly when a delocalized 
electron is added to the neutral molecule. This is due 
to the fact that there is a near degeneracy between un- 
occupied onefold and twofold states. The spin-orbit in- 
teraction leads to a strong mixing between these states 
which, because of the orbital components, enhances the 
anisotropy. 

In the following we discuss quantum transport 
through a Mni2 molecule weakly coupled to metal- 
lic electrodes. Electron tunneling between leads 
and the molecule is described by the Hamiltonian 
■^T = Ea,/ I]fe,p*aittLcfc(p) + H.C., where c\ (ck) cre- 
ates (destroys) an electron in orbital k in the molecule, 
and ta depends on the tunneling barrier between lead 
a = L, R (left or right) and the molecule. The leads, 
at electrochemical potential /ic, are described by the 
independent-electron Hamiltonian Ha = 
where a^^j (^a/) creates (destroys) a quasiparticle of 
quantum number /. We also take into account the 
work function for the external leads. As a result, the 
charge state populated at zero gate and bias voltages 
is the neutral one, and not the anion as it would 
seem from Table I. The sequential tunneling current 
is obtained from a master equation for the occupation 
probabilities of the molecule many-body states. The 
transition rate between two many-body states via 
tunneling of an electron from lead a into the molecule, 
is proportional to /(i?(p', fc-, TVg/) — E{p,ki,NQ) — 
Ma) X |(p';fc;,/c^,...,fc5v^j4(p)|p;A:i,fc2,...,fcjvQ)|2, 
where f{E) is the Fermi distribution function and 
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FIG. 1. (color online), (a) Lowest lying spin multiplets for the 
three charge states of the molecule; the small energy difference 
between states M and — M is due to transverse anisotropy 
terms coming from 4th-order spin-orbit contributions and our 
choice of the 25* -I- 1 values of {9m , (jJM ) • See text after Eq. ((2]| . 
(b) Differential conductance of a Mni2 SET as a function of 
bias and gate voltage, for a symmetric double junction, with 
gate capacitance equal to 1/20 of the total capacitance, at zero 
temperature. Numbers -1,0 and -1-1 denote the excess charge 
on the molecule and label regions of Coulomb blockade. In 
regions (—1,0), (0,-1-1) and ( — 1,0,-1-1) transport is possible 
via transitions between different charge states. 



FIG. 2. (color online), (a) Differential conductance as a func- 
tion of bias at Vg — —20 V. Parameters are as in Fig. [T] 
Solid (dotted, offset for clarity) lines: calculation based on 
SDFT (giant-spin model). Inset of (a) shows a zooming of 
the onset of SDFT transport due to transitions between the 
ground-state (GS) spin multiplets of the neutral and anionic 
molecule. Visible peaks in the main plot correspond to tran- 
sitions between the GS spin multiplet of the neutral and the 
first three excited multiplets of the anion, (b) Spin multiplets 
involved in the transport. The transitions between the states 
generating the dominant peak in (a) are indicated by arrows 
joining the states. 



i?(p, fcijiVg) is the energy of state |p; fci, fc2, • ■ • , fcjVg), 
modified by the bias Vf, = (^l — fJ,R)/e and gate voltage 
Vg. We then solve numerically the master equation in 
steady-state and obtain the current as a function of Vb 
and Vg. Figure [TJb) shows the differential conductance 
G — dl/dVb as function of Vb and Vg. The calculations 
are done at zero temperature. We choose equal coupling 
of the molecule to the two leads; the gate capacitance 
is equal to 1/20 of the total capacitance of the system. 
Three CB stability diamonds are visible, corresponding 
to the three different charge states Q — —1,0, +1, where 
transport is blocked. In region indicated by (—1,0) 
[respectively (0, -|-1)], current flows through transitions 
between anionic (cationic) and neutral states. In region 
(—1,0,-1-1) all three charge states are present. The 
additional lines parallel to the GS-GS transitions are 
due to transitions between excited states. In Fig. [TJb) 
we can also see two lines, indicated by arrows, that 
correspond to a decrease in the current with increasing 
Vb (NDC). These lines give the bias at which, for a given 
Vg, anionic states become occupied in the (0,-1-1) region. 
NDC in Mni2-SET has been observed experimentally 
i. 

For a better understanding of transport just above the 
CB gap, in Fig. [Ifa) we plot the differential conductance 
as a function of Vb, for Vg — —20V. Transport in this 
region is due to transitions between the spin multiplets 
of the neutral and anionic molecule. The conductance 
peak spectrum displays a rich fine structure, with peak 
spacing on the order of 0.1-1 meV, which corresponds 
to that seen in experiment 0, Q ■ The first set of peaks 



at the very onset of transport is caused by transitions 
between the GS spin multiplets. Surprisingly, the con- 
ductance in this region is very small, G < 10~*e^//i, as 
shown in the inset; it is practically invisible for transi- 
tions between the low- lying states (large \AI\ and |M'|) 
and slightly larger for transitions between higher-lying 
states (small \M\ and |Af'|). As we argue below, this is 
caused by the very small overlap of the orbital parts of 
the many-body wave-functions of the two GS spin mul- 
tiplets. The second set of peaks in Fig. [21[a), Vb > 1.485 
V, corresponds to transitions between the GS spin mul- 
tiplet of the neutral molecule and the first three excited 
spin multiplets of the anion. This cluster of resonances is 
largely determined by the first excited spin multiplet of 
the anion, since reaching this multiplet opens up trans- 
port also via other multiplets. In particular, the domi- 
nant peak seen in the figure is due to transitions between 
the lowest-energy states of the GS spin multiplet of the 
neutral molecule and the first excited spin multiplet of 
the anion, as shown in Fig. [2jb). 

In order to shed light on the interplay between orbital 
and spin-selection rules, we compare the SDFT-based 
calculation with the giant-spin model. Within this spin 
model, transitions are possible only between states whose 
spin differs by 1/2 (spin-selection rule), with transition 
rates given by Clebsch-Gordan (CG) coefficients [24*1. In 
the computation of the conductance, we include the GS 
spin multiplet of the neutral molecule, and the GS and 
first three excited spinmultiplets of the anion. The con- 
ductance for the giant-spin model is shown with the dot- 
ted line in Fig. ^a). The first 11 peaks in the inset cor- 
respond to subsequent transitions between states M = 



4 




FIG. 3. (color online). Matrix elements for transitions 
between anionic and neutral charge states, (a) Giant-spin 
model; (b) SDFT results for transitions from neutral GS spin 
multiplet to anionic first excited multiplet; (c) SDFT results 
for transitions from neutral GS to anionic GS spin multiplets. 

±S, ±{S -1),...,0 and M' = ±S', ±{S' - 1), ... , 1/2, 
where S' — S + 1/2. The intensity of these peaks de- 
creases monotonically with decreasing \M\, which is dif- 
ferent from the SDFT-based conductance. The more 
complicated set of peaks at Vb > 1.485 V in Fig. ^a.) 
resembles the analogous cluster of peaks for SDFT and 
has the same interpretation. We examine the matrix el- 
ements giving the neutral-to-anion transition rates. Fig- 
ure[3]Ja) shows the results for the giant-spin model, where 
the matrix elements are proportional to the CG coeffi- 
cients and the spin selection rule |M — M'\ — 1/2 is 
strictly obeyed. Figure EJb) shows that the same spin- 
selection rule is approximately satisfied by the SDFT 
matrix elements of the transitions between the GS spin 
multiplet of the neutral molecule and the first excited 
spin multiplet of the anion. In particular, the orbital 
part of the wave functions does not modify substantially 
this condition. In contrast. Fig. [3] (c) shows that the 
SDFT matrix elements for the transition between the two 
GS spin multiplets are different: the effect of the spin- 
selection rules is now overridden by space selection rules, 
which suppress most of the transition rates near |M| = 5 
and lead to a vanishing conductance. Furthermore, the 
GS matrix elements close to the diagonal behave dif- 
ferently as a function of |M| for the two models: they 
decrease with \M\ for SDFT and increase for the spin 
model, which is reflected in the conductance [Fig. ^a)] 
for Vb < 1.485 V. Based on Figs. [2Ia-c), we expect the 
giant-spin model to agree better with the SDFT calcula- 
tion for transitions involving the first excited-state spin 
multiplet of the anion. Indeed, Fig. [UJa) shows that for 
bias voltages Vb > 1.485 V the two models yield qual- 
itatively the same conductance. The small matrix ele- 
ments in Fig. ^c) are also the cause of the NDC seen in 
Fig. [2b) along the line separating the transport regions 
(0,-1-1) and (—1,0,-1-1): when Vb becomes large enough 
to access the anion GS multiplet, the system remains 
trapped in these states due to their small connection to 
the neutral states. Thus, the current decreases. 



In conclusion, we presented a microscopic study of the 
GS properties and low-energy spin excitations for the 
neutral and charged Mni2 molecular magnet, based on 
SDFT. Resonances in the tunneling conductance are gov- 
erned both by spin and spatial selection rules. The latter 
ones play a key role in determining the relative contribu- 
tion to transport of various spin multiplets, and can lead 
to NDC. The orbital properties of the spin states pro- 
vided by SDFT are essential to build a correct effective 
spin model and interpret the transport experiments. 
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Note added - After this manuscript was submitted, an- 
other Letter (2^ related to our work appeared. 
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